Enhanced production of direct photons in Au + Au collisions at v /s~^~=200 GeV 

and implications for the initial temperature 
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The production of e + e" pairs for m e + e - < 0.3 GeV/c 2 and 1 < Pt < 5 GeV/c is measured 
in p + p and Au + Au collisions at y/s NN = 200 GeV. Enhanced yield above hadronic sources is 
observed. Treating the excess as photon internal conversions, the invariant yield of direct photons 
is deduced. In central Au + Au collisions, the excess of direct photon yield over p + p is exponential 
in transverse momentum, with inverse slope T — 221 ± I9 stat ± 19 syst MeV. Hydrodynamical models 
with initial temperatures ranging from Ti n i t —300-600 MeV at times of —0.6-0. 15 fm/c after the 
collision are in qualitative agreement with the data. Lattice QCD predicts a phase transition to 
quark gluon plasma at —170 MeV. 

PACS numbers: 13.85.Qk,25.75.Cj,12.38.Mh,21.65.Qr 



Experimental results from the Relativistic Heavy Ion 
Collider (RHIC) have established the formation of dense 
partonic matter in Au + Au collisions at ^/snn = 200 
GeV The large energy loss of light quarks and glu- 
ons [2| as well as that of heavy quarks 3] indicates that 
the matter is very dense. The strong elliptic flow of 
light 0, [fl and charmed 0] hadrons indicates rapid ther- 
malization. Such a hot, dense medium should emit ther- 
mal radiation [||; the partonic phase is predicted to be 
the dominant source of direct photons with 1 < pt < 3 
GeV/c in Au + Au collisions at RHIC Q- 

Observation of thermal photons will allow determina- 
tion of the initial temperature of the matter. However, 
the measurement of direct photons for 1 < pr < 3 GeV/c 
is notoriously difficult due to a large background from 
hadronic decay photons. Direct photons contribute only 
~ 10% above the background photon yield Q ■ In general, 
any source of high energy photons can also emit virtual 
photons, which convert to low mass e + e _ pairs. For ex- 
ample, gluon Compton scattering (q + g — > q + 7) has 
an associated process that produces low mass e + e~ pairs 
through internal conversion (q + g —* q + "f* — ► q + e + e~). 
Consequently, we search for "quasi-real" virtual photons, 
which appear as low invariant mass e + e~ pairs. 

The relation between photon production and the asso- 
ciated e + e~ pair production can be written as 



2a 1 

37T m e . 



Ami 
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Here a is the fine structure constant, m e and m ee are the 
masses of the electron and the e + e~ pair respectively, and 
S is a process dependent factor that goes to 1 as m ee — > 
or m ee <C pt- Equation |T]) also describes the relation 
between the photons from hadron decays (e.g. ir°,r] —> 
77, and u> — > 77T°) and the e + e~ pairs from Dalitz decays 
( 7r°,?7 — > e + e _ 7 and to — ► e + e~ir°). For 7r° and 77, the 



factor S is given by S = |F(m 2 J| 2 (l - r -^f [10], where 

Mh is the meson mass and F{m\ e ) is the form factor. 

The factor S for a hadron h is zero for m ee > Mh. 
We exploit this cut-off to separate the direct photon sig- 
nal from the hadronic background. Since 80% of the 
hadronic photons are from tt° decays, the signal to back- 
ground (S/B) ratio for the direct photon signal improves 
by a factor of five for m ee > 71^0 = 135 MeV/c 2 , thereby 
allowing a direct photon signal that is 10% of the back- 
ground to be observed as a 50% excess of e + e _ pairs. 

In this Letter we present the analysis of e + e~ pairs for 
m ee < 0.3 GeV/c 2 and for 1 < p T < 5 GeV/c in Au + Au 
and p + p collisions at y/s NN = 200 GeV recorded during 



2004 and 2005, respectively. The PHENIX detector [11 1 
measures electrons in the two central arms, each cover- 
ing I At; I < 0.35 in pseudorapidity and 7r/2 in azimuthal 
angle. The Au + Au analysis 0, [l2| uses 8xl0 8 mini- 
mum bias (Min. Bias) events corresponding to 92.2j^|' q% 
of the inelastic Au + Au cross section. The beam-beam 
counters and zero degree calorimeters provide the Min. 
Bias trigger, as well as the centrality selection 13J. The 
p+p analysis [3] uses 43 nb _1 of data recorded using the 
Min. Bias trigger and 2.25 pb _1 of single electron trig- 
gered data. Helium bags in both runs reduced the total 
conversion material, including the beam pipe, to ~0.4% 
of a radiation length. 

All electrons and positrons with px > 0.2 GeV/c are 
combined into pairs. Pairs from photon conversions in 
the detector material are removed by a cut on the orienta- 
tion of the pair in the magnetic field [§] . The combinato- 
rial background is computed by mixing events and is sub- 
tracted [E El. The S/B ratio is -0.2 (at m ee =0.3 GeV) 
to -1.5 (at m ee = 0.1 GeV/c 2 ) for p T > 2 GeV/c and 
0.05 to 0.4 for 1 < px < 2 GeV/c. There are two sources 
of correlated background: two e + e~ pairs from a meson 
decay and correlated hadrons decaying into two e + e~ 
pairs, either within the same jet or in back-to-back jets. 
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FIG. 1: (color online) The measured e + e~ pair invariant mass 
distributions. The px ranges are shown in the legend. The 
solid curves represent an estimate of hadronic sources; the 
dashed curves represent the uncertainty in the estimate. 



The magnitude of the correlated background, about 10% 
of the signal in p + p, is determined from the like-sign 
pair data and subtracted after correcting for acceptance 
differences between like and unlike-sign pairs [14j. We 
correct for electron reconstruction efficiency, and in p + p 
for trigger efficiency, determined as a function of mass 
and pair px using a GEANT-based Monte Carlo simula- 
tion of the PHENIX detector. 

Figure [T] shows the mass spectra of e + e~ pairs in 
p + p and Au + Au collisions for different ranges of 
pair pt, comparing to a "cocktail" of hadron decays 
calculated using a Monte Carlo hadron decay generator 
based on meson production measured by PHENIX [9(. 
Detector resolution is included in the cocktail calcula- 
tion. The op en charm contribution, calculated with 
PYTHIA [16], is also included but is negligible in this 
kinematic range. The cocktail is normalized to the 
data for m ee < 0.03 GeV/c 2 ; the absolute normaliza- 
tion agrees with the data within a 20% systematic un- 
certainty [IT 



14] 



The "knee" beginning at m ee ~ 
0.1 GeV/c 2 corresponds to the ir° cut-off, leading to an 
80% reduction of background above this point. The p+p 
data are consistent with the background for m ee > M^o 
at lower px, but reveal a small excess over the back- 
ground at higher px- A much greater excess is observed 
in Au + Au indicating enhanced production of virtual 
photons. 

Internal conversion of direct photons is a possible 
source of the excess. Little contribution from other 
sources of e + e~ pairs is expected in this mass region since 
7r + 7r~ — ► e + e+ can only contribute for m ee > 2M V . Al- 
though PHENIX has observed a strong enhancement of 
e+e~ pairs for 0.15 < m ee < 0.75 GeV/c 2 in Au + Au, 
it peaks at low px and decreases rapidly with increas- 
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FIG. 2: (color online) Electron pair mass distribution for 
Au + Au (Min. Bias) events for 1.0 < pr < 1.5 GeV/c. 
The two-component fit is explained in the text. The fit range 
is 0.12 < m ee < 0.3 GeV/c 2 . The dashed (black) curve at 
greater m ee shows /(m ee ) outside of the fit range. 



ing pr@] with a different mass distribution than that 
observed at high px- 

Figure [2] shows that the mass spectrum for m ee < 
0.5 GeV/c 2 and px > 1 GeV/c is well described by 
the cocktail plus internal conversion photons. The flat 
mass spectrum of the excess above the cocktail at this 
px shows no significant indication of low-mass enhance- 
ment 9j. Thus, we treat the excess entirely as internal 
conversion of photons and deduce the real direct photon 
yield from e + e _ pairs using Eq. (fTj). 

We fit a two-component function f(m ee ) — (1 — 
r)f c (m ee )+r / d i r (m ee ), to the mass distribution. / c (m ee ) 
is the shape of the cocktail mass distribution (shown in 
Fig- H]), /dir(m-ee) is the expected shape of the direct 
photon internal conversion, and r is the fit parameter. 
We assume that the form factor for direct photons is 
F( m ee) — 1> as one would expect from a purely point-like 
process. For direct photons from parton fragmentation 
or from hadronic gas, F[vr? ee ) may be greater than one. 
If we arbitrarily set the form factor in /dir(w ee ) to be the 
same as that in f^{m ee ), r would decrease by ~ 10%. 

For each px bin, f(m ee ) is fit to the data for m\ ow < 
m ee < 0.3 GeV/c 2 with m low = 0.08,0.1,0.12 GeV/c 2 ; 
r is the only fit parameter. Figure [2] shows fdwijnee) 
and f c (m ee ) together with a fit result for Au + Au (Min. 
Bias) data for 1.0 < px < 1.5 GeV/c. For higher px bins, 
\ 2 /NDF is near 1.0; fits to centrality separated data also 
give good x 2 '/NDF. 

Therefore, we focus on the uncertainties that can cause 
distortions in the mass distribution, namely (i) the par- 
ticle composition in the hadronic background, (ii) the 
background (from mixed events and correlated pairs), 
(iii) the geometric acceptance due to detector active ar- 
eas, and (iv) the efficiency corrections. These were stud- 
ied by Monte Carlo simulation. The mass spectrum is 
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FIG. 3: (color online) The fraction of the direct photon com- 
ponent as a function of pr- The error bars and the error band 
represent the statistical and systematic uncertainties, respec- 
tively. The curves are from a NLO pQCD calculation (see 
text). 



distorted within the systematic uncertainties, and the 
fitting procedure is applied to the distorted spectrum to 
determine the systematic uncertainties in r. The sys- 
tematic uncertainty due to the variation of m\ ow is also 
included. The dominant uncertainty is the particle com- 
position in the hadronic cocktail, namely the 77/ 7r° ratio 
which is 0.48 ±0.03(0.08) at highp T for p + p (Au + Au) 
based on PHENIX measurements 17] . This corresponds 
to a ~ 7% (~ 17%) uncertainty in the p + p (Au + Au) 
cocktail for 0.1 < m ee < 0.3 GeV/c 2 . Other sources 
cause only a few percent uncertainty in the data to cock- 
tail ratio. 



Figure [3] shows the fraction r of the direct photon com- 
ponent determined by the two-component fit in (a) p + p 
and (b) Au + Au (Min. Bias). The curves represent 
the expectations from a next-to-leading-order perturba- 
tive QCD (NLO pQCD) calculation pi]. For p + p, 
the curves show the ratio da NLO 
da^ LO { PT 

NLO pQCD calculation and da™ cl (pT) is the inclusive 



7 (p T )/da l ™ l (p T ), where 
is the direct photon cross section from the 



photon cross section. For Au -I- Au, the curves represent 
^(p T )/dN™ x {p T ), where T AA is the Glauber 



T AA da^ LO < 



nuclear overlap function and dN!^ cl (pT) is the inclusive 
photon yield. The three curves correspond, from top to 
bottom, to the theory scale fi = 0.5 pt, Pt, and 2 px, 
respectively, showing the scale dependence of the theory. 
While the fraction r is consistent with the NLO pQCD 
calculation 18] in p + p, it is larger than the calculation 
in Au + Au for p T < 3.5 GeV/c. 

The direct photon fraction r in Fig. [3] is converted to 
the direct photon yield as dN dir (p T ) = r x dA incl (p T ). 
The inclusive photon yield gL/V 1i1c1 (pt) for each pt bin 
is determined from the yield of e + e~ pairs for m ee < 
0.03 GeV/c 2 using Eq. ([1]). Here we use the fact that in 



AuAu Min. Bias x10 
AuAuO-20%x10 2 
AuAu 20-40% x10 
P+P 
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FIG. 4: (color online) Invariant cross section (p + p) and in- 
variant yield (Au + Au) of direct photons as a function of pr- 
The filled points are from this analysis and open points are 
from [3, |2C| . The three curves on the p + p data represent 
NLO pQCD calculations, and the dashed curves show a modi- 
fied power-law fit to the p+p data, scaled by Taa- The dashed 
(black) curves are exponential plus the Taa scaled p + p fit. 
The dotted (red) curve near the 0-20% centrality data is a 
theory calculation |7J. 



this mass range the process dependent factor S is unity 
within a few percent for any photon source. 

Figure [4] compares the direct photon spectra with pre- 
viously measured direct photon data from [l^ . |2(]| and 
NLO pQCD calculations [ijj]. The systematic uncer- 
tainty of the inclusive photon (14% from the uncertainty 
in the e + e~ pair acceptance correction [l^]) is added in 
quadrature with the systematic uncertainties of these 
data. The p + p data are shown as an invariant cross 
section using da — a'^dN. 

In this analysis we have converted the yield of excess 
e + e+ pairs to that of real direct photons using Eq. ([1]), as- 
suming S = 1. This implies = ^-^-dn~ l . Thus the 

1 am ee 6tt m ee i 

yield of the excess e + e~ pairs for 0.1 < m ee < 0.3 GeV/c 2 
before the conversion can be obtained by multiplying the 
direct photon yield by a factor of|^log|22 = 1.7xl0 -3 . 

The pQCD calculation is consistent with the p+p data 
within the theoretical uncertainties for px > 2 GeV/c. A 
similarly good agreement is observed for tt° [2l]]. The 
p + p data can be well described by a modified power-law 
function ( J 4 p?) (l-|-py/&) _n ) as shown by the dashed curve 
in Fig.UJ The Au + Au data are above the p+p fit curve 
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scaled by Taa for px < 2.5 GeV/c, indicating that the 
direct photon yield in the low px range increases faster 
than the binary NN collision scaled p ± p cross section. 



TABLE I: Summary of the fits. The first and second errors 
are statistical and systematic, respectively. 



centrality 


dN/dy(p T > lGeV/c) 


T(MeV) 


X 2 /DOF 


0-20% 


1.50 ±0.23 ±0.35 


221 ± 19 ± 19 


4.7/4 


20-40% 


0.65 ±0.08 ±0.15 


217 ± 18 ± 16 


5.0/3 


Min. Bias 


0.49 ±0.05 ±0.11 


233 ± 14 ± 19 


3.2/4 



We fit an exponential plus the TAA-scaled p ± p fit 
function {Ae- pT / T + T AA x A pp (l + p T /b)- n ) to the 
Au ± Au data. The only free parameters in the fit are 
A and the inverse slope T of the exponential term. The 
systematic uncertainties in T are estimated by chang- 
ing the p ± p fit component and the Au + Au data 
points within the systematic uncertainties. The results 
of the fits are summarized in Table HI where A is con- 
verted to dN/dy for pt > lGeV/c. For central collisions 
T = 221 ± 19 stat ± 19 syst MeV. Using, instead, a power- 
law function (oc p^ n ) to fit the p ± p spectrum yields 
n = 5.40 ±0.15, and T Au Au = 240 ±21 MeV. If the direct 
photons in Au ± Au collisions are of thermal origin, the 
inverse slope T is related to the initial temperature Ti n i t 
of the dense matter. In hydrodynamical models, Tj n i t is 
1.5 to 3 times T due to the space-time evolution [22l |. 
Several hydrodynamical models can reproduce the cen- 
tral Au ± Au data within a factor of two [9]. These as- 
sume formation of a hot system with initial temperature 
ranging from Ti n ; t = 300 MeV at thermalization time 
t = 0.6 fm/c to T init = 600 MeV at r = 0.15 fm/c [13]. 
As an example, the dotted (red) curve in Fig. 2] shows a 
thermal photon spectrum in central Au ± Au collisions 
calculated with T init = 370 MeV Q. 

In conclusion, we have measured e + e~ pairs with 
m ee < 300 MeV/c 2 and 1 < p T < 5 GeV/c in p + p 
and Au ± Au collisions. The p + p data show a small 
excess over the hadronic background while the Au + Au 
data show a much larger excess. By treating the excess 
as internal conversion of direct photons, the direct pho- 
ton yield is deduced. The yield is consistent with a NLO 
pQCD calculation in p + p. In central Au ± Au colli- 
sions the shape of the direct photon spectrum above the 



TAA-scaled p ± p spectrum is exponential in px, with an 
inverse slope T = 221 ± 19 stat ± 19 syst MeV. Hydrody- 
namical models with T init ~300-600 MeV at r ~0.6-0.15 
fm/c are in qualitative agreement with the data. Lattice 
QCD predicts a phase transition from hadronic phase to 
quark gluon plasma at ~170 MeVQ. 
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